Entanglement conditions and polynomial identities 
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We develop a rather general approach to entanglement characterization based on convexity proper- 
ties and polynomial identities. This approach is applied to obtain simple and efficient entanglement 
conditions which work equally well in both discrete as well as continuous-variable environments. 
Examples of violations of our conditions are presented. 
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One of the most fascinating features of Quantum Me- 
chanics is the phenomenon of entanglement. Despite 
its long history, originated with the seminal Einstein, 
Podolsky and Rosen paper [l[ in 1935, many questions 
about this phenomenon arc still open. The basic prob- 
lem is the ability to distinguish entangled states from 
non-entangled (i.e. separable) ones. The first complete 
characterization of entangled states has been given in Q 
in the form of a hierarchy of conditions whose first level 
is the famous Peres- Horodccki positivity of partial trans- 
position condition [3, 0] (another approach has been de- 
veloped in Q). This hierarchy is a numerical procedure, 
with all restrictions and limitations of a numerical algo- 
rithm in a setting where an analytical description of the 
problem at hand is desirable. That is why the existence 
of a complete numerical characterization does not elimi- 
nate the need in analytical conditions. 

A real-valued function f(g), defined on the set of all 
density operators, is called convex if the inequality 



(1) 



is valid for all states g n and for all numbers p n with 
the normalization condition ^2 n p n = 1- The simplest 
example of a convex function is the average value (A) g 
for a fixed operator A, for which the condition (JXJ) be- 
comes equality. If F : R — > R is a convex function of a 
real argument then F((A)g) is also convex. For example, 
el and (A) 2 ~ are convex functions of g (for fixed A). 
In fact, {A) 1 - is convex if n = 2n' is an even number. 

A real- valued function g(g) is called concave if the in- 
equality, opposite to the one given by Eq. ([1]), is valid 



(2) 



In other words, a function is concave if it increases on 
mixtures of states. Examples of a concave function are 
the average (A)§ and the variation cr^(g) = (A 2 ) — (A) 2 

of the observable A. Note, that the average (A)g, being a 
linear function of the quantum state, is both convex and 



concave. A trivial special case is the operator A = a ■ 1 
of multiplication by a constant a, or simply a constant 
operator, for which (A)§ = a is both convex and concave. 

After these preliminary definitions we can formulate 
the following lemma, which plays a central role in our 
work: 

Lemma 1. Let f(g) be a convex function and g{g) be 
concave. If the inequality 



f(0n) < g(g n ) 



(3) 
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is valid for some states g n , n = 1,2,..., then it is also 
valid for all their mixtures ^ n p n gn, where p n ^ and 
T, n Pn = 1- 

Proof. The proof is rather simple: for a given mixture 
g = ^2 n p n gn we have the inequalities 

f(e) < J2pnf(g n ) «S ^p n g{g n ) «S g{g)- (4) 



Roughly speaking, the left-hand side does not increase 
on mixtures, but the right-hand side does not decrease, 
which completes the proof. □ 

As a consequence of this lemma, we obtain the follow- 
ing result. If the inequality is valid for all factorizablc 
states Q n = gn® Q b n ® ■ ■ ■ ® Qn then it is also valid for 
all separable states. A violation of this inequality is a 
sufficient condition for entanglement. This characteri- 
zation is complete in the sense that for any entangled 
quantum state g one can find such a pair of functions / 
and g which satisfy Lemma [1] (with g n being separable 
states) and for which f(g) > g(g). In fact, as it has been 
proven in for any entangled state g there is an ob- 
servable E such that (E)+ ^ for all separable states f, 
but for which (E) s < 0. Let us take f(g) = (0 • i) s = 
and g(g) = (E)§, then / is convex (and concave) and g 
is concave (and convex). According to the choice of E 
the inequality (0) is satisfied for all separable states, but 
f(o) > g{o) f° r the given entangled state g. 

Though this characterization is general and complete, 
it is not very useful in practice since it is not clear how 
to choose appropriate functions / and g. Thereby, the 
main difficulty to get efficient entanglement conditions 
using Lemma Q] is to find pairs of functions / and g, 
one of which is convex and the other is concave, and 



2 



which satisfy the inequality ([3]) for all factorizable states. 
Moreover, it must be possible to violate this inequality. 
Finding such pairs of function is the main result of our 
work. 

Before we formulate our first result, we prove a cou- 
ple of simple lemmas. The following one deals with the 
variances of factorizable states. 

Lemma 2. Factorizable quantum states g = g a <E> g b sat- 
isfy the inequality o~ AB (<? a ® Q b ) ^ o~A{g a )o~B{Q b ) for arbi- 
trary observables A and B of different degrees of freedom. 

Proof. In fact, for the factorizable state g we have the 
relation (AB) g = (A) ga (B) ^ and analogously for the 
squares of observables, so we have 



o-ab(q) - (A 2 )^(B 2 )g b {A)UB)\ 



(5) 



Using the identities of the form (A 2 ) = a A (g) + (A) 2 , we 
obtain the relation a AB (g) = a 2 A {g a )a 2 B (g b ) + . . ., where 
the dots on the right-hand side of this equality stand for 
the sum a A (g a ) (B) 2 gb + (A) 2 - a a B (g b ). Since this expres- 
sion is always non-negative we get the desired inequal- 
ity. □ 

The next lemma concerns quadratic polynomials, non- 
negative on the right half of the real line and puts some 
condition on the coefficients of such polynomials. 

Lemma 3. If the inequality a A 2 — bX + c 0, is valid 
for all X ^ 0, where a,b,c^ 0, then b 2 ^ 4ac. 

Proof. In fact, if a = then b must be zero too, in other 
case this inequality, which in this case reads as — bX + c 
0, is violated for sufficiently large A. If a > then the left- 
hand side of the inequality can be written as aX 2 — bX+c = 
a(X - b/{2a)) 2 +c- 6 2 /(4a), and for A = b/(2a) > it 
must be non-negative, from which we obtain the desired 
inequality b 2 ^ 4ac. □ 

Now we are ready to prove one of our main results. 

Theorem 1. For arbitrary observables A, A' and B, B' 

of different degrees of freedom the inequality 



o-abo-a'b> > ^\([A,A'][B,B'])\ 
is valid for all bipartite separable quantum states. 



(6) 



Proof. Using Lemma [2j we can estimate the quantity 
°~ab (s) + °~\'B' (s) f° r a factorizable state g = g a <g> g b as 
follows: o 2 AB + a 2 A , B , ^ o 2 A (g a )a 2 B (g b ) + a 2 A ,(g a )a 2 B ,(g b ). 



The right-hand side of this inequality is bounded from 
below by the product 2a A (g a )a A >(g a )a B (g b )a B >(g b ). Ap- 
plying the Heiscnbcrg uncertainty relation for the ob- 
servables A and A' to the product of the first two vari- 
ances (expressed as a A (g a )a A , (g a ) ^ (1/2)|([1, A'])\§ a ) 
and the same relation for the observables B and B' to the 
product of the last two variances and taking into account 



the identity {[A, A']) s « ({B, B']) gb 
finally get the inequality 



([A,A'][B,B'}) § , we 



'AB 



i 



o-a>b>(q)>^\([A,A>}[B,B'}} s \ 



(7) 



for factorizable states. As we have already mentioned, 
the variance o~ 2 AB {g) is a concave function of g (for fixed 
observables A and B) and the quantity | (O) g | is convex, 
so we can apply Lemma [T] and conclude that the inequal- 
ity ([7]) is valid for all bipartite separable states. Now lets 
scale the observables A and B by a non-negative factor 
\/A, A ^ 0: A \/Ai, B \fXB. After scaling we get 
that any bipartite separable quantum state satisfies the 
inequality o\ B X 2 - (l/2)|([i, A'] [B, B']) \X + a\, B , ^ 
for all non-negative A. Applying Lemma 02 we get the 
desired inequality ([6]). □ 

The proof of this theorem shows that the inequality 
given by Eq. ([6]) immediately follows from Lemma [2] for 
a factorizable state. But one cannot use Lemma [2] for 
a general separable state, since the left-hand side of the 
inequality ©, which is a product of two convex func- 
tions, may not be convex, as the product x 2 (x — l) 2 of 
two functions of real argument illustrates. Both factors 
(parabolas) are convex but their product is not. That is 
why we need the trick with the scaling of observables. 

It is interesting to compare the inequality ([6]) and the 
Heisenberg uncertainty relation for the observables AB 
and A'B'. The latter gives us the following lower bound 
for the product of the dispersions of the observables un- 
der study: oaboa'B' > (l/2)|([AB,i'B'])|. Note that 
this inequality is valid for all bipartite quantum states. 
Our inequality ([6]) gives another lower bound for the 
same product of dispersions, but this bound is guaran- 
teed to be valid only for separable states (there can be 
entangled states which also satisfy this bound). If, for a 
given state and observables, we have |([^4, > 
2\([AB 7 A'B'})\, then the inequality ([6]) gives a stronger 
restriction then the Heisenberg relation so it may be pos- 
sible to violate it (this means to calculate the product of 
dispersions and compare it with (1/4)| ([A, A'} [B, B']) |). 
In other case, the Heisenberg uncertainty relation pro- 
hibits violations of the inequality ©. As we will see 
below, our inequality can be violated in many important 
cases, both discrete as well as continuous variable ones. 
As a measure of violation we use the quantity 



V 



\([A,A'][B,B'] 
4er abo A'B' 



(8) 



According to the inequality ([6]), V < 1 for all separable 
states, so V > 1 is a sufficient condition for entanglement. 

Our first example of violation of the inequality (J6j) is 
for the especially simple choice of the observables A = x ai 
A' = p a , B = p and B' = Xb, the position and momen- 
tum operators. In this case the inequality ([6]) reads as 



0~xz>0~ri 



1 



(9) 
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We have \([AB,A'B'})\ = | Im(a 2 -S 2 )|, where a and b are 
annihilation operators of the first and the second mode 
respectively, and in many cases this quantity is zero. Let 
us take the state 



+00 

E 

n=0 



c n |2n, 2n) 



(10) 



with real coefficients c„, where \n) are the Fock states. 
One can easily check that for this state we have (x a Pb) = 
(ftaXb) = 0. To compute the dispersions note that 
(xlfb) = (pl%b) and that 

x 2 Jl = \{{2a^a+l){2b%+l)-a% 2 -aH^)+. . . , (11) 

where the dots stand for the terms that do not contribute 
to (x^pl). The product of the dispersions reads as 



Vxp&px = ^ + Q(co, Ci, . . 



(12) 



where Q = Q(cq,Cx, . . .) is a quadratic form of the coef- 
ficients c„, n = 0, 1, . . . 

+00 

Q = (M2n + 1)<$ - (« + l)(2n + l)c„c„ +1 ) . (13) 

n=0 

The matrix of this quadratic form is given by 



C = 



( 1/2 

1/2 6 3 

3 20 15/2 

15/2 42 14 

14 72 

v 



(14) 



Let us consider a truncated (N + 1) x (N + 1) matrix 
Cjv, which corresponds to the truncated state (|10|) (where 
Crj = if n > N). The minimal value of the expression 
c t Cnc for c = (cq, ci, . . . , cat) with the normalization 
condition |c| = 1 is the minimal eigenvalue A m j n of the 
matrix C/y and the corresponding eigenvector gives the 
coefficients Cfe, k = 0,1, . . . , N. With Mathematics we 



have obtained that A m ; n tends to the value 



-0.04495 



as n — > +00. This corresponds to the maximal value 



1/4 



(1/4) + A r 



1.2192. 



(15) 



For the truncated state \4>2) = co|00) + Ci|22), where Cq 
and Ci are real and c\ + c\ 



1, we have a xp 



px 



(1/4) + 6c 2 — cqC\. If we take the coefficients such that 
< ci < co/6 then we get a xp a px < 1/4. The violation 
V as a function of Co (assuming that c\ = yl — Cq) takes 
the maximum for Co ~ 0.997, ci w 0.077 and the maximal 
value is V ~ 1.197. This means that even the simplest 
truncated state ^2) gives a rather good violation. 

Now consider a mixture of the state (1101) with vacuum 



(16) 



£=2#><VI + (1-P)|00)(00|, 



where ^ p ^ 1. It is easy to check that for this mixed 
state we have 



PxpVpx = ^ +PQ{co, Cl, ■ • ■)• 



(17) 



For p = 0, when the state (j!6p is simply vacuum, there 
is no violation, as it must be. But for any < p ^ 1 we 
have 



Vmaxb) 



1/4 



1/4 + pA n 



>1, 



(18) 



since A m ; n < 0. We see that even a tiny fraction of the 
entangled state (|10p in the mixture (TIB")) can be detected 
with the condition 

Let us choose two orthogonal quantum states |0) and 
|1), (0|1) = 0. Define the Pauli-like operators s x , s y , s z 
and sq as follows: 



|0)(1| + |1)(0|, S v = -i(|0)<l|-|l)(0|), 
|0)(0|-|1)(1|, « = |0)<0| + |1)(1|. 



(19) 



If the system is two-dimensional, then |0) and |1) form a 
basis in the state space of the system and Sq is the identity 
operator §0 = 1. The operators f| 1Q[) satisfy the well- 
known spin-like commutation relations [s x ,s y ] = 2is z , 
[s y ,s z ] = 2is x , [s z ,s x ] = 2is y , the equalities s x =s y — 
J 2 = so and the anticommutator relations {s x ,s y } — 
{s x , s z } = {s y , s z } = 0, where {A, A'} = AA' + A' A. 

We show that our condition can detect entanglement 
of any pure bipartite state of the form 



If) 



*\lp )\ip ) + fllpl)^ 



(20) 



where \tpo) and \ip±) is a pair of orthogonal states of the 
other mode of the system and |a| 2 + |/3| 2 = 1. Since ([20]) 
is, in fact, the Schimdt decomposition of the state |^), 
this state is separable if and only if either a = or j3 = 0. 

One can verify that (s x s b x ) = -{s^s b y ) = 2Re(a/3*) 
and (s a J b y ) = (s a y s b x ) = -2Im(a/3*), (sjs^) = (sgsg) = 1. 
Now consider the operators 



A = s° cos + Sy sin 9, A' = s% cos 6' + s° cos 9' , 
B = s x cos 77 + s y sin 77, B' = s b x cos 77' + s b y sin 77'. 

Then wc have 

(AB) = 2Re(a/3*e l(e+,)) ), 
(A'B') = 2Rc(a/3*e l(e '+"' ) ), 



(21) 



(22) 



[A, A'} = 2isin(9' - 9)§ a z and [B,B'} = 2?;sin(77' - rj)s b z . 
Let us take the angles 9 = — arga, 77 = arg/3, 9' = 9 + 
7r/2, 77' = ?/+7t/2, and we have (AB) = -(A'B') = 2\a/3\, 
[A, A'] = 2is% and [B,B'\ = 2is b z . Then the inequality 
© reads as 1 — 4|a/3| 2 ^ 1, which is obviously violated 
unless a = or {3 = 0, that is unless the state ([20| is 
factorizable. Note that in this way one can verify en- 
tanglement of any pure state with the Schmidt number 
larger than 1. 



4 



The inequality ([6]) can also be violated by the squeezed 
vacuum state defined as 



+00 



\S) = VT^J2x n \n,n) 



(23) 



n=0 



with real parameter A, |A| < 1. If for both modes we take 
the operators A = J2n=o \ 2n )( 2n + M + \ 2n + l )( 2n \ and 
B = En=o -i(\ 2 n)(2n + 1\ + \2n + l)(2n|), we get 



V = 



1 + A S 
1-A 5 



(24) 



We see that in this case we have V > 1 for A 7^ and 
V —> +00 when A — > 1. 

Let us consider the inequality obtained in Q , which in 
the bipartite case reads as 



(25) 



(Ab - A'B') 2 + (A& + A'B) 2 

«S ((A 2 +A' 2 )(B 2 +B' 2 )). 

It has been proven in Q that this inequality cannot be 
violated in the case of quadratures. For two level systems 
and dichotomic observables this inequality becomes 



(AB ~ A'B') 2 + (A& + A'B) 2 < 4, 



(26) 



which has been proven to be valid for all bipartite two- 
level quantum states In it has been shown that 
the multipartite analogue of the inequality (|25p can be 
violated for a large number of modes (more than 10). The 
inequality ([25]) easily follows from our approach with the 
help of the identity 

(ab - a'b') 2 + (ab' + a'b) 2 = (a 2 + a' 2 )(b 2 + b' 2 ), (27) 

which expresses the multiplicativity of the norm of the 
complex numbers (consider the product of two complex 
numbers a + ia! and b + ib'). Now for any factorizable 
state we have 

(Ab - A'b') 2 + (Ab' + A'b) 2 = (Ab) 2 + (Ab') 2 
+ (A'b) 2 + (A'b') 2 < ((A 2 + A' 2 )(b 2 + B' 2 )). 

(28) 

As it has already been noted, the left-hand side of this 
inequality is convex and the right-hand side is concave 
and Lemma Q] shows that the inequality (|2"5"|) is valid for 
all separable states. The multipartite analogue of the 
inequality (f25j) can be obtained using the identity for 
product of more then two complex numbers. Complex 
numbers is not the only numerical system whose norm 
is multiplicative. Inequalities based on the multiplicativ- 
ity of the norms of quaternions and octonions have been 
obtained in Q. 

Nothing prevents us from using the identity (|27|) in the 
other direction and derive the inequality 



((AB - A'B') 2 ) + ((A& + A'B) 2 ) 

> (Ab) 2 + (Ab') 2 + (A'b) 2 + (A'b') 2 , 



(29) 



which after some algebraic manipulations can be trans- 
formed to the following form: 

vab + vab'+oa'b+oa'B' > \([A,A'][B,B'})\. (30) 

It is clear that this inequality follows from ©. 

The theorem Q] can be straightforwardly generalized to 
the multipartite case. 

Theorem 2. For arbitrary observables A\, A\, . . . , A n , 

A' n of different degrees of freedom the inequality 

<TA,...A n <T A{ ... A ' n >^\([A 1 ,A' 1 ]...[A n ,A' n })\ (31) 

is valid for all completely separable multipartite quantum 
states. 



It is easy to see that the Bell state 



1 



|s„) = - 7 =(|o...o) + 11...1)) 



(32) 



violates this inequality for A^ = si k \ A' k = which 
in this case reads as 



3 y ' 



0-x...xVy... V > \(s\ 



(33) 



If n is even then the state \B n ) is an eigenstate of the 

operators s^ . . . s x n , s^ . . . and si 1 ^ . . . si™^ with 
eigenvalues ±1, so the left-hand side of the inequality 
(f33|) is zero and the right-hand side is 1. 

Our next result is based on the Ramajuan identities. 

Theorem 3. For any observables A, A' and B, B' of 

different degrees of freedom the inequality 

(Ab + Ab' + A'B) n + (Ab' + A'b + A'B') n 

+ (Ab - A'B') n ((Ab' - A'B) n )+ (34) 

((A'b + A'b' + Ab) 71 ) + ((A'b' + Ab + Ab') 71 ) 

with n = 2 and n = 4 is valid for all bipartite separable 
states. 

The proof is based on the Ramanujan identity 

(ab + ab' + a'b) n + (ab' + a'b + a'b') n 

+ {ab-a'b') n = (a'b + a'b' + ab) n (35) 
+ (a'b' + ab + ab') n + (ab' - a'b) n , 

which is valid for n = 2, 4, see [l(| , and proceeds as above. 
In this case we can also get two inequalities, but one can 
easily see that the two are, in fact, the same, after an 
appropriate permutation of the variables. 

To illustrate the possibility to violate the inequalities 



3j), let us take A = s a x , A' = s£, B = s b x , B' = s b y and 



s b B' = sb 

apply these inequalities to the simplest Bell state 

1 



|B 2 ) = -==(|00> + |11)). 



(36) 
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The inequality (|34|) for n — 2 reads as 



\ s x S x • s x s y "r s y s xl • \ s x s y ' s y s x < °y S yl 

sbx2 <8-6(^), 



and for n = 4 as 



(S^ + ^ + ^4> 4 + (s a J b y + s«S b x + s a y s b ) 4 
+ (S a J b x -'s a y s b y ) 4 ^34-32(s a J b z ). 



(37) 



(38) 



It is easy to check that (s%s b ) = ~{s^s b ) = (s a z s b z ) = 1, 



(s%s b ) = (s^s b x ) = 0. For the inequality ([37]) we have 
6 = l 2 + (— l) 2 + 2 2 < 2, which is clearly wrong. The 
inequality (jSHJ) becomes 18 = l 4 + (-1) 4 + 2 4 < 2, which 
is also wrong, so both these inequalities can be violated. 

In conclusion, we have developed a general approach to 
entanglement characterization and obtained simple and 
powerful entanglement conditions, which can be used for 
discrete as well as for continuous-variable systems. Ex- 
amples of violations of these conditions are presented. 
Multipartite generalization of some of our results are also 
discussed. 
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